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Abstract: Let 9 > 0. We consider a one-dimensional fractional Ornstein-Uhlenbeck process defined 
as dXt = —0 Xtdt + dBt, t > 0, where _B is a fractional Brownian motion of Hurst parameter 
H £ (i, 1). We are interested in the problem of estimating the unknown parameter 9. For that pur- 
pose, we dispose of a discretized trajectory, observed at n equidistant times ti — zA„, i — 0, . . . ,n, 
and Tn — nA„ denotes the length of the 'observation window'. We assume that A„ — > and T„ — > oo 
as n — > oo. As an estimator of 9 we choose the least squares estimator (LSE) 6'„. The consistency 
of this estimator is established. Explicit bounds for the Kolmogorov distance, in the case when 
H e |), in the central limit theorem for the LSE 0„ are obtained. These results hold without 
any kind of ergodicity on the process X . 

Key words: fractional Ornstein-Uhlenbeck processes, discrete-time observation, least squares esti- 
mator, Kolmogorov distance, central limit theorem, Malliavin calculus. 

1 Introduction 

In this paper we consider a fractional Ornstein-Uhlenbeck process X = {Xt,t > 0). That is, it solves 
the linear stochastic differential equation 

Xq = xo; dXt ^ ~9 Xtdt + dBt, t>0, (1.1) 

where xq Q R, B = {Bt, i > 0) is a fractional Brownian motion with Hurst parameter H G (5,1) and 
> is an unknown parameter. 

Assume that the process X is observed equidistantly in time with the step size A„: ti = iA„,i = 
0, . . . , n, and T„ = nA„ denotes the length of the 'observation window'. The purpose of this paper 
is to study the least squares estimator (LSE) 9n of 9 based on the sampling data At . , z = 0, . . . , n. 
The LSE 0„ is obtained as following: 9n minimizes (formally) 



2 

n I I 



i=l 

where ti = iA„, i — 0, . . . ,n,. Thus 9n is given by 

_ Yli=iXu-i{Xt^ - At,_ J 

^nl^i^l^ti.i 

Also, by using (jl.lL we arrive to the following formula: 



where 

U^ = Xt^ - Xt^_, + 9AnXt^_, , i^l,...,n. 



(1.2) 



(1.3) 
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The parametric estimation problems for fractional diffusion processes based on continuous-time 
observations have been studied e.g. in [51 [T31 [T2] via maximum likelihood method. Recently, the 
parametric estimation of continuously observed fractional Ornstein-Uhlenbeck process defined in 
(jl.ip is studied in [HIH], in the case when H G (i, 1), by using the least squares estimators. 

From a practical point of view, in parametric inference, it is more realistic and interesting to 
consider asymptotic estimation for fractional diffusion processes based on discrete observations. 
There exists a rich literature on the parameter estimation problem for diffusion processes driven by 
Brownian motions based on discrete observations, we refer to and [12] . In this paper, we focus 
our discussion on the LSE and the fractional Ornstein-Uhlenbeck process case. 

In general, the study of the asymptotic distribution of any estimator is not very useful for practical 
purposes unless the rate of convergence of its distribution is known. The rate of convergence of the 
distribution of LSE for some diffusion processes driven by Brownian motions based on discrete time 
data was studied e.g. in [7]. To the best of our knowledge there is no study of this problem for the 
distribution of the LSE of the unknown drift parameter in equation p.ip . Our goal in the present 
paper is to investigate the consistency and the rate of convergence to normality of the LSE 9n defined 
in (fra . 

Recall that, iiY,Z are two real- valued random variables, then the Kolmogorov distance between 
the law of Y and the law of Z is given by 

dKoi(r, Z) = sup \P{Y <z)-P{Z<z)\. 

— oo<2;<oo 

Let us now describe the results we prove in this work. In Theorem 13.31 we show that the consistency 
of On as A„ — )■ and 7t.A„ — oo holds true if H G (^j !)■ When H e |) we use the Malliavin 
calculus, the so-called Stein's method on Wiener chaos introduced by [S] and the technical Lemmas 
13.41 and 13. SI proved respectively by [6, and to derive Berry-Esseen-type bounds in the Kolmogorov 
distance for the LSE 9n (Theorems 13.61 and 13.71) . 

We proceed as follows. In Section 2 we give the basic tools of Malliavin calculus for the fractional 
Brownian motion needed throughout the paper. Section 3 contains our main results, concerning the 
consistency and the rate of convergence of 0n- 

2 Preliminaries 

In this section we describe some basic facts on the stochastic calculus with respect to a fractional 
Brownian motion. For more complete presentation on the subject, see [9| and [1]. 
The fractional Brownian motion {Bt, t > 0) with Hurst parameter HE (0, 1), is defined as a centered 
Gaussian process starting from zero with covariance 

Rnit, s) = E{B,B,) = \ {t'" + s'" \t s\'") . 

We assume that B is defined on a complete probability space {fl, T , P) such that J- is the sigma-field 
generated by B. By Kolmogorov's continuity criterion and the fact 

EiBt-Bsf = \s-t\^"; s, t> 0, 

we deduce that B has Holder continuous paths of any order j < H. 

Fix a time interval [0, T]. We denote by H the canonical Hilbert space associated to the fractional 
Brownian motion B. That is, H is the closure of the linear span £ generated by the indicator functions 
l[o,t], t G [0;^] with respect to the scalar product 

(l[o,t]j l[o,s]) = Rnit, s). 
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The application if ^ £ — > B((p) is an isonietry from £ to the Gaussian space generated by B and it 
can be extended to H. 

U H S (5,1) the elements of H may be not functions but distributions of negative order (see pD]). 
Therefore, it is of interest to know significant subspaces of functions contained in it. 
Let \H\ be the set of measurable functions (p on [0, T] such that 



dudv < 00. 



M^^^:^H{2H-l) f f \^iu)Mv)\\u-v\^ 
Jo Jo 

Note that, if eln], 

E{B{ip)B{^)) = H{2H ~1) [ I ip{u)il){v)\u~v\'^"-^dudv. 

Jo Jo 

It follows actually from [10] that the space \'H\ is a Banach space for the norm || .|| |^| and it is included 
in In fact, 

i2([0, T])cLTr ([0, T])c\H\cH. (2.1) 

Let C^(K",M) be the class of infinitely differentiable functions / : M" — > R such that / and all 
its partial derivatives are bounded. We denote by S the class of smooth cylindrical random variables 
F of the form 

f{Biip,),...,B{cp„)), (2.2) 
where n > 1, / € Cg°(R",]R) and ipi,...,ipn e H. 

The derivative operator D of a smooth cylindrical random variable F of the form (|2.2|) is defined as 
the "H- valued random variable 

In this way the derivative DF is an element of L^{il;H). We denote by D^'^ the closure of S with 
respect to the norm defined by 

\\F\\l2^Ei\\Fr) + E{\\DF\\l). 

The divergence operator S is the adjoint of the derivative operator D. Concretely, a random variable 
u £ L^{ft;7i) belongs to the domain of the divergence operator Domd if 

E\{DF,u)H\<Cu\\F\\mn) 
for every F £ S. In this case 6{u) is given by the duality relationship 

E{F5{u)) = E{DF,u)^ 
for any F £ D^'^. We will make use of the notation 



6{u) 



r 

I UsSBs, u e Dom{S). 
Jo 



In particular, for h £ H, B{h) = 5{h) = hsSBg. 

Assume that H £ (i, 1). If u e D^'^dHD, u belongs to DomS and we have (see [SJ Page 292]) 



^( ll^w|lrHi®i«i 
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where the constant ch depends only on H. 

As a consequence, applying (|2.ip wc obtain that 



E{\S{u)n < CH \\EiuW 1 



E \\Du\ 



L77([0,T]2) 



(2.3) 



For every n > 1, let Hn be the nth Wiener chaos of B, that is, the closed linear subspace 
of L^{V,) generated by the random variables {Hn{B{h)),h e %, \\h\\-H = 1} where is the nth 
Hermite polynomial. The mapping /„(/i'^") = n\Hn{B{h)) provides a linear isometry between the 



symmetric tensor product H®" (equipped with the modified norm ||.||?^0" = :^^|MIh®") ^^'^ ^n- 
For every f,gG H®" the following multiplication formula holds 

£^(/n(/)/„(ff))=n!(/,5)««„. 

Finally, It is well-known that L^(i7) can be decomposed into the infinite orthogonal sum of the 
spaces Tin- That is, any square integrable random variable F E L^(i7) admits the following chaotic 
expansion 

oo 

F^E{F) + Y,In{fn), 
n=l 

where the /„ S H®" are uniquely determined by F. 

We will make use of the following theorem proved in [5]. 

Theorem 2.1 (Nourdin-Peccati). Let F e D^''^ has zero mean and N ^ A/'(0, 1). Then 



dKoi{F,N)<JE {l-<DF,-DL-^F>u) 



(2.4) 



Moreover, if F = Iq{f) with q > 2 and f £ H®'?, then < DF, -DL-^F ^WDFW^ and therefi 

dKoi{F,N) < 



\ 



E 



^--WDFWli 



(2.5) 



Fix T > 0. Let f,g : [0,T] — > R are Holder continuous functions of orders a S (0,1) and 
/? € (0, 1) with a + (3 > 1. Young [14] proved that the Riemann-Stieltjes integral (so-called Young 
integral) fsdgs exists. Moreover, ii a — (3 £ 1) and : — > K is a function of class C^, the 

integrals ^{fu, gii)dfu and ^{fu, gu)dgu exist in the Young sense and the following change of 
variables formula holds: 



Hft,9t) ^ <l>{fo,go) + Tr;{fu,gu)dfu+ Tp (/«,5«)'^5«, < t < T. 



df 



dg 



(2.6) 



As a consequence, a H e (i,l) and (ut, t e [0,r]) be a process with Holder paths of order 

a e (1 — if, 1), the integral /^^ UsdBs is well-defined as Young integral. Suppose moreover that for 
any t e [0,r], ut G and 



P 



T 

Jo 



\DsUt\\t- s\'^"^'^dsdt < oo 



Then, by [T], m € DomS and for every t e [0, T], 

[ UsdBs = [ UsSBs+H{2H -1) [ [ D^Url 
Jo Jo Jo Jo 



.\2H-2 



drds. 



(2.7) 
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In particular, when ip is sl non-random Holder continuous function of order a G (1 — i/, 1), we obtain 

rT pT 



f ifisdB, ^ f ^sSBs = B{ip). (2., 
Jo Jo 



In addition, for all i/i S 

E\l ipsdBs i ^psdBs] = H{2H -1) [ [ ip{u)ip{v)\u - v\'^"~'^dudv. (2.9) 



/ ^sdBs [ ^PsdBs] =H[2H-l) [ [ 
^v/o Jo J Jo Jo 



3 Asymptotic behavior of the least squares estimator 

Throughout this paper we assume H G (|, 1) and 9 > 0. Let us consider the equation (|l.ip driven 
by a fractional Brownian motion B with Hurst parameter H and 9 is the unknown parameter to be 
estimated for discretely observed X. The linear equation (II. ip has the following explicit solution: 

-St I ^ I / „6ls. 



Xt = e-'* [xo + J e'^'dB.j , t > 0, (3.1) 

where the integral can be understood either in the Young sense, or in the Skorohod sense, see indeed 
Let us introduce the following process 

:= / e^'dBs, t > 0. 







We shall use the notation a„ s3 &„ to means that there exists a positive constant c(xo, 9, H) (depend- 
ing only on xo,9 and H) such that, 

sup |a„|/|6„| < c{xq,9,H) < oo. 

n>l 

We define the following sequence, which will be used throughout, 

a„ = H{2H - 1) / / e-''^"u^"~^dudt, n > 0. (3.2) 



JO Jo 

We shall be using the following lemmas several times. 

Lemma 3.1. Let H G (j, 1), let 9 > 0, and let Q!„ be the sequence defined by iS.'t^} . Then 



JO 



XsdB,+an, n>0, (3.3) 



and 



In particular, an ^ r„. 



lim = e^-^"HT{2H). (3.4) 
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Proof. By (|2.7p . we have 



( XsdBs = [ XsSBs + H{2H -1) [ [ DsXt\t - s\'^"-^dsdt 
Jo Jo Jo Jo 

= [ " X,SBs + H{2H -1) [ " [ e-^^*-'\t^ sf^-^dsdt 
Jo Jo Jo 

= [ " XJBs + H{2H ~l) f " f e-^^u^^-^dudt 
Jo Jo Jo 



Jo Jo 

XsSBs + a„ 







On the other hand, 



an ,. H(2H-l) r e„ off 2 , , 

hm — = hm — ^ / / e-^^'u^^'^dudt 



Tn 







/>oo 

= H{2H - 1) / e-''"f|2H-2^^ 
= e^-^"HT{2H). 

Thus, the proof is finished. □ 

Lemma 3.2. Assume H G (1/2, 1) and 9 > 0. Then, there exists a constant c > 0, depending only 
on xq, 9 and H , such that 

supE (Xf) < c<oo. 
t>o 

Proof. For any t > 0, we have 

E{^f) = H{2H-1) f f e'^e'y\x-yf"-^dxdy 
Jo Jo 

= 2H{2H-l)f dye^y r dxe^'^iy - x) 
Jo Jo 

= 2H{2H-\) I dye^'^y C dze 
Jo Jo 

= H{2H ^1) J^ze^'^'z^"'^ J' dye 

= H{2H-1) I dze-^'^z 
Jo 



2H-2 



t „20t „2ez 



HT{2H)2et 



Q2H 

By combining p.ip and p.5p . we obtain that for any i > 



e''\ (3.5) 



This proves the claim. □ 



3.1 Consistency of the LSE 

The next statement provides consistency of the LSE 0„ of 6. 

Theorem 3.3. Assume H E (1/2, 1) and 9 > 0. Then, if An — > and nA„ — oo as n oo, we 

have 



Proof. From (|1.3p , we can write 



in probability as n ^ oo. 



(3.6) 



Let < p < 1. We have 



Pi 



P 



< P 



9 " 



i=l 

■= + j2{n). 

We begin studying the term ji{n). We have, 



>p(l-p) 1 +P 



9Ar. 



> p 



P 



P 



rv.. ' ^ 



> p(i - p) 



,=1 Ju^i 



+- 



"■n Jo 



XtSBt 



>P(1-P) 



< P\ — 



-P 



> 3/5(1 -P) 



-p 



i=l 

XtSBt 




> 3^(1 -P) 



ji,i(»^) + Ji,2(n-) + ji,3(»^)- 
For the term ji^i(n), by using Lemma 13.21 and the fact that for every t > 

Xu_,-Xt = (e-«*--e-''*)(xo+6,-J+e-^*(6,-, -6) 

we obtain 

n n 

Y,E\Xt^_,[U.-{Bu-Bu_,)\\ = E^ 



(3.7) 



i=l 



i=l 



Xu_, I 9{Xt^_, - Xt)dt 

ti-i 
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< 



■1=1 -'^.-1 

Making the change of variables s = t — we obtain 



^ 1/2 

i=i 

^/ (3.8) 

i=l 

i=i 

+E/ ^-'"{E{[i,-^,^_,r)y/'dt 

J2 [1 - e-^(*-*-^)]rft 
i=i 

„■_. ■Jl.i-^ 



E / [1 - e-^(*-**-^)]di = n / [1 - e-^^rfs 

,_i Jti-i Jo 



,J^^"[l-e-'^]ds 



n 



where the last estimate comes from the fact that "[1 — e ^'^Jds/A^ — > 0^/2 as A„ — >• (thanks to 

t-fi-1 ' t-ti-1 



I'Hopital rule). On the other hand, by the change of variables u = ^ , v = ^ and s = t—ti-i 



E/ e-^*(^^([6-6.-J^))'^'rfi 

= ^/H{2H-l)f2 r dte-"' ( f dye'y f dxe^''\x - y\''"-^] 

= ^H{2H-l)Y, r rfte-^(*-**-) i f dye'^y-'^-^^ f dxe'^--'^-^^\x - y\^"-A 

= ^yH{2H-l)J2 dt{t - ii_i)^e-^(*-*-^) (^^' dt;e^(*-*-i)" due'^'-**-^^^\u - v\^"-^^ 



1/2 



^ n I s^ds 







Hence, we obtain that 



i=l 
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which leads to 
Consequently, 

= p(— 



J- n. 



A" 



For the term ji,2{n), from p.7p . we have 



E 



Using the inequality (I2.3p . we can write 



i ^^^^ 



< CH 



CH 



CH 



^0 



-I 



E(^ 



1=1 



dsdt 



H 



JO 



E 



i=l 



H 



/.A 



(1-e 



-e(t-t,_i) 







= chtl 



H 



(3.9) 



(3.10) 



because " [l - e'^"] dv/A;r' ^ 9^/" /{I + 1/H) as A„ ^ (thanks to I'Hopital rule). 
Similarly, we have 



E 



i=l 



= E 



< CH 



i=l 



"'0 



1/H 



H 



dsdt 



9 



H 



CH 



CH 



CH 



ChU 



1=1 
T cT ^ 



H 



E 



l[t._i,t](s)l[t._i,til(Oc'sc'i 



/ c?f / c?se 



-e{t-s)/H 



h 



JO 



H 



" I I dv due-'-/" 

J<3 



= CHu" y 



[1 - e-^"/^] 

e/H 



H 



dv 



Thus, 



E 



Therefore, 

As consequence, 



*=1 -^i-i 



t,-i 



A? 



For the term Ji,3(?t.), by settmg 



Ft 



iTn Jo 



XtdBt 



we have 



ji.sH = ^ — 

\ Ctn 



< 9r 



Tn 



XtSBt 

I 



> 3^(1 - P) 



since from [3], 



E{F}J ^ A{9,H), asn^oo. 



(3.11) 



(3.12) 



(3.13) 



(3.14) 



(3.15) 
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where 



11-4H Tj2 



H\AH - 1) 



T{2H) 



2 r(2i?)r(3-4i7)r(4i/- 1) 



r(2 - 2H) 



Finally, by combining (|3.10l) . p.l2p and (|3.14p . we conclude that 



1 



1 



(3.16) 



W - pW 

Consequently, to achieve the proof of Theorem 13.31 it remains to estimate the term j2 [n) . We have 



P 



= P 



a. 



■E^t.-i >p 

-A„x:^L.-i 



i=l 



> P 



j2,i{n) + j2,2{n). 









> p/2 j + P 1^ 


— r" x^dt - 1 


>pI^ 





We first estimate J2,l(?^)• Similar argument applied in p. 81) yields 



< 



< 



■E 



20sup,>o(i^[^t'])^/'^ 



1/2 



A2 + nA^+i 



which implies that 



sg A, 



p 



H 



i=l 



>p/2 



A^ 



We now study j2,2{n). Applying the change of variable formula p.6p leads to 

XtdBt = -xl + 2e / Xt^dt. 

Jo 

Combining p.3p and p.lSp we obtain 



(3.17) 



(3.18) 



261 



X2dt-2a„ = 2 / XtSBt-X^ + xl- 



Hence 



J^n Jo 



X^dt - 1 



•i^ri Jo 



Xt5Bt-—{X^ -xl). 



(3.19) 
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We deduce from ((3?T9|) and ((3?T5|) together with that 



P 



< P 



^ P 



x?dt - 1 



>p/2 



(3.20) 



$1 



1 

Oin JO 

1 

P^Tn pT„ 

1 



> p/8 + P 



>p/2 



2a„ 



>p/8 



+ 



pTn 



1 



Therefore, from (I3.17P and p.2ip . we obtain 



A" 1 

J2{n) ^ ^ + 



(3.21) 

(3.22) 
□ 



P ' P^Tn 

Finally, combining p.l6p and p. 221) . the proof of Theorem 13.31 is done. 
3.2 Rate of convergence of the LSE 

This paragraph is devoted to derive Berry-Esseen-type bounds in the Kolmogorov distance for the 
LSE 9,1 of 6. We first recall the following technical lemmas. 

Lemma 3.4 ([BJ Page 78]). Let f and g be two real-valued random variables with g ^ P-a.s. Then, 
for any 6 > 0, 

dKoi(^^,N^ <dj,,,{f,N) + P{\g-l\>S) + S. 

where TV- A/'(0, 1). 

Lemma 3.5 ([H Page 280]). Let Y and Z be two real-valued random variables. Then, for any rj > 0, 



dKoi {Y + Z,N)< dKoi {Z, N) + P{\Y\ >ri) + 



/27r 



where N ^ N{0,1). 
Define 



An 



OTnjEiF^ 



n > 1, 



where Q!„ and Ft^ are defined by (j3.2l) and (I3.13p . respectively. 
The following result provides explicit bounds for the Kolmogorov distance, in the case when H £ (i, |), 
between the law of A„-\/2^(0„ — 9) and the standard normal law 
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Theorem 3.6. Let {6,r]) e (0, 1)^. If H e (i, |) then, for some constant c > depending uniquely 
on Xq, and H, we have: for any n > 1, 



d-Kol 

< c 



XnVTnien -0),N 



A" 



rj ' ' '5 nA„(5^ 

where iV-7V(0, 1). 

Proof. Fix [S^rj) G (0, 1)^. From (|1.3p and Lemma we obtain that 



< du 



XnVTni0n-9),N 
1 



Y,Xu-iU^,n\ +P 



:= Ji(7i) + J'iin) + (5. 
We first study the term Ji{n). Using Lemma [3751 we obtain 



>(5 +(5 



Ji(n) 



1 



= dy. 



1 



T„JE{F^ ) .=1 



< P 



/T^^/E{F^J 



i=l 



> 



P 



+dj. 



/T^jE{Fl^ 



{Xt,_, - Xt,)SBt 



> 



Xt5Bt,N 



/T^^E{F^JJo 
:= Ji,i(n) + Ji,2("-) + -^i,3("-) + 
By using dSl]), (pUj) and (151^ we deduce that 



2n 



r,VT^^E{Fl) 



(3.23) 
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To achieve the estimation of Ji{n) its remains to estimate Ji,3(n). We have 



1 



XtSBt = 



1 



Thus, by using Theorem [Q and the fact that E (||£'-Pr„ ||^) = 2E {F^J, we obtain 



T^JE{F^ ) -Jo 



XtSBt, N\ = dK 



-.,N 



< 



\ 









2 \ 2- 


E 








[h 













2E{F}J 



Moreover, from [3], we have 



E 



{WDFtJI - E\\DFtJI) 



S3 t: 



E 



IDFtJI^EWDFtJIY 



Thus 



Consequently, by combining p.23p and p.24p we have 

7? 



(3.24) 



(3.25) 



FinaUy, via the same arguments as in the estimation of j2 [n) , we obtain the following upper bound 
of J2{n): 



A" 1 

J2{n) ^ + -— ^, 



which completes the proof of Theorem 



(3.26) 
□ 



Let us apply now Theorem 13.61 to the particular case rf = V^nA^i^^ and 5 ~ A" , where < a < 
(3 < H which ensure that 7iA,\+2" — ^ oo and nA}^^^ — >■ as n — ?> oo. 



Theorem 3.7. Let < a < l3 < H such that nA^+2" -!• oo and nA]+'^P as n oo. If 

^2' 4) 



H G (i, |) then, for some constant c > depending uniquely on xq, 9 and H , we have: for any 



n> 1, 



dKoi {XnVT^idn ~d),N) < c (^Af-^ + {nA,,r"-^ + ^ nAi^^^ + A^-^ + + A^ 

In particular, as ?i —> 00 



/l^(0„-0)^AA(O,a|,) 



wh 



ere (t| = {AH - 1)6 [l 



r(3-4g)r(4g-i) 
r(2-2H)r(2ff) 
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